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ABSTRACT:  The  non-linear  steady-state  heat  conduction  equation 

which  arises  in  the  theory  of  thermal  explosions,  was  solved 
by  Frank-Kamenetzky  for  the  case  of  the  semi-inf inite  slab, 
and  by  Chambre  for  the  solid  cylinder  and  sphere  geometries. 
Solutions  for  the  case  of  the  :o!low  cylinder  are  presented 
and  it  is  shown  tnat  from  these,  the  slab  and  solid  cylinder 
solutions  can  be  deduced  as  special  cases.  The  design  of 
large  rocket  grains,  where  self-heating  may  result  in 
"spontaneous"  ignition  during  manufacture  or  storage.  Is 
considered. 
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These  computations  were  performed  under  Task  N0L-S6-2d-2-l-56, 
"Fundamentals  of  Solid  Propellant  Ignition  and  Burning". 

While  present  rocket  grain  webs  are  restricted  to  thicknesses 
which  render  spontaneous  combustion  of  conventional  propellants 
unlikely,  the  possibility  must  not  be  ignored.  In  particular, 
propellants  with  a  low  activation  energy,  low  thermal  con¬ 
ductivity,  high  density  and  high  rate  of  decomposition  are 
suspect.  Within  the  assumptions  made,  the  results  of  this 
study  are  believed  to  be  accurate  and  are  the  responsibility 
of  the  author  and  the  Fuels  and  Propellants  Division  of  the 
Naval  Ordnance  Laboratory. 


W.  W.  WILBOURNE 
Commander 
Captain,  USN 


A.  E.  ABLARD 
^oy  direction 


il 


NAVORD  Report  4267 


TABLE  OR  CONTENTS 

Page 

INTRODUCTION .  1 

MATHEMATICAL  MODEL . 1 

SEMI— INR'I NITE  SLAB  LIMIT . 5 

SOLID  CYLINDER  LIMIT . 7 

ENGINEERING  APPLICATIONS .  8 

REFERENCES . 8 


111 


NAVORD  Report  426? 


THE  SOLUTION  OF  THE  STEADY  -STATE  HEAT  CONDUCTION 
EQUATION  WITH  CHEMICAL  REACTION  FOR  THE  HOLLOW  CYLINDER 

INTRODUCTION 


Let  a  combustible  undergo  an  exothermic  chemical  reaction 
with  heat  loss  to  the  walls  of  the  containing  vessel  where  this 
loss  is  assumed  to  take  place  via  a  conductive  process  inside 
the  combustible  volume.  The  criterion  for  thermal  explosion 
is  that  the  energy  liberated  during  the  reaction  must  be 
greater  than  that  lost  through  the  surfaces.  The  condition 
when  the  heat  loss  exactly  compensates  energy  production  in 
the  combustible  is  described  by  the  steady  temperature  state. 


MATHEMATICAL  MODEL 

The  steady  state  heat  conduction  equation  can  be  written 
as 

XV*T  =  -/>ClZ.x|>(-£)  ,  u) 

where  T  is  the  combustible  temperature,  X  the  thermal  con¬ 
ductivity,  fi  the  density,  GL  the  heat  of  reaction,  Z.  the 
frequency  factor,  E  the  activation  energy,  R  the  gas  con¬ 
stant,  and  V2  the  Laplacian  operator.  It  is  assumed  that  the 
reaction  is  unimolecular .  The  solution  of  Eq .  (1)  for  the  semi- 
infinite  slab  immersed  in  an  isothermal  bath  at  temperature * 

,  under  the  assumption*  T  -T^ «  was  first  given  by 
Frank-Kamenetzky  ( D .  With  the  same  assumption,  Chambr£(2) 
showed  that  the  solution  of  Eq.  (1)  for  solid  cylindrical  and 
spherical  geometries  can  be  obtained  in  terms  of  known  functions. 

Consider  a  hollow  combustible  cylinder  of  inner  radius 
a  and  outer  radius  b  whose  isothermal  surfaces  are  at 
temperatures  "TJ  and  TJ,  respectively.  To  solve  Eq.  (1)  assume 
chat  T-Tk  «Tfe  .  Then  to  first  ord^r  in  T-T^  > 

*■  The  unpublished  numerical  integration  of  Eq.  (1)  for  the 
slab  and  solid  cylinder  performed  at  the  U.  S.  Naval  Ordnance 
Laboratory  has  shown  that  this  assumption  is  very  good  for  the 
range  E/rt  £  fc  . 


1 


NAVORD  Report  42  67 


r  rk  1  +  0.  Tb  ) 

£  .  K  v  f. 

Setting  &  =  )  •)  2-=-g-  , 

Equation  (1)  becomes  in  cylindrical  coordinates, 

0  +  iii  =  _fexp^ 

where  f*  is  the  radial  space  coordinate  and 


(2) 


(3) 


(4) 


Let  Ism  whenrsa  •  Equation  (3)  is  to  be  solved  subject  to  the 
boundary  conditions 


at  the  outer  surface 


(5) 


(6) 


0=0,  l»l  - , 

and  at  the  inner  surface, 

S’  =  ,  *  =  m  <  I 

Noting  the  method  of  solution  by  Chambre,  the  variables 

yi  =  2Z  exp  (6)  ;  w  =  a  >  (7) 

are  introduced  and  Equation  (3)  becomes 

duj  _  _  S*  (8) 

d  7  2,+  uj 

wnose  general  integral  is 

(9) 


u>z*-4u>+D  =  -  2.  $  Yi 


The  constant  D  (which  was  zero  in  Chambre' s  work)  is  in  general 
different  from  zero.  Rewriting  Equation  (9)  in  terms  of  the 
original  variables,  there  results 
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which  when  compared  with  Equation  (3),  yields 


The  solution  of  this  equation  is 

■0*  =  -  2  Xn  (Ae'  +  ^  +•  Bz'  *  ) 

where  ^  :  l“—  . 

4 

Boundary  condition  (5'  gives 


(10) 


(11) 


(12) 


B  =  I- A 


so  that 


e  =  -z  [a*1*1  +  G-a)z'~*] 

Boundary  condition  (6)  is  used  to  evaluate  A  ,  yielding 


A=  ifm 


-  IV\ 
.  .-2* 


I  -  m 

where  5  X«K*  Equation  (12)  still  contains  the 

parameter  •R  whose  relationship  to  S“  will  be  found.  Sub¬ 
stituting  Equation  (14)  into  Equation  (3)  gives 

S’  =  8  &Z  A(l-A)  .  (it 

Eliminating  A  between  Equations  (15)  and  (16)  yields 

f  -  g  %z  O-rm^'Kvm*"1  -ha2*) 
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which  is  the  desired  relationship  between  £  and  “fe  for  given 
m  and  Y  «  From  Equation  (17)  it  can  be  shown  that  Tor  given 
m  and  t  ,  J  has  a  maximum  critical  value,  Sc  ,  for  which  the 
steady-state  solution  is  still  possible.  When  the  value  of  S 
as  defined  by  Equation  (U)  exceeds  7C  ,  then  thermal  explosion 
will  occur. 


Introduce  the  new  parameter 

oC  =  vn  *  (18) 

so  that  Equation  (17)  becomes 

f  =  ivi*  <!-«*)*  (19) 

Letting  O  »  it  is  found  that  the  root,  <*c  ,  of  the 

equation  * 

+*  _  -  jl _ +  .jy - 

<x4-  I  oc  \  imf*/  Y w-m  m<x-jr  ^2°) 

makes  S  a  maximum,  Sc  ,  for  given  m  and  Y  .  From  Equation 
(18)  is  obtained  the  critical  value  of  ft  ,  •  Using 

Equations  (llj.).  (15)  *  and  (18),  the  critical  temperature 
distribution,  -0-c  ,  is 

%  -  Z (M,1) -  i V [(l-l «t) \l‘ ^].(a  1 

The  critical  temperature  gradient  is  given  by 

£ 

ae£  _  _  z  {  +  (22) 

2  (w»-y#c,)*z  c  +(y-  ro<xc)*c 

At  the  boundaries. 
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If  -f  and  'ft  represent  the  energy  flux  and  energy  transfe] 
rate  per  unit  axial  length,  then 


<47* 


fr-X^  ;  i  =  -in \ril  ■ 

dr  or 


(25) 


Equation  (25)  <*  n  be  written  as 

Fr-J?  ;H=-2^, 

dz  d  Z 


where 


;H  = 


(26) 


(27) 


XRT*^  '  "  irrXRV 

The  position,  i0  ,  of  the  maximum  critical  temperature,  ^  ^ 
which  occurs  where  -  O  t  Is  * 

I 


K  = 


(*-mrwxc)«« 

rv\~y<x. 


2.X. 


|  I  »**""• 

Substituting  Equation  (28)  into  Equation  (21)  gives 


(28) 


A-  -  _ i  f  m-Y0tg  ,  ^Cy-bxxc-)qfc  ,  l-Bc 


(29) 


Consider  the  case  when  the  inner  and  outer  surfaces  are 
at  the  same  temperature,  i.e.,  =  l  •  Table  1  gives  the 

values  of  <Xc,yt,  Ft.m  *  F»>»  »  W|Jh,  ,Hls,  ,  2.  ,  and 
»  for  various  values  of  m  . 

C.  **i  JiX 


SEMI -IN UNITE  SLAB  LIMIT 

Prom  Table  1  it  is  seen  that  as  I  ,  o(c  is  finite  and 
»  oo  •  Equations  (23),  (24),  (26)  and  (29)  give 
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Jk«v  =  ^  - 

»*i  c 


-|  ;  Jkvw  &  -  J2*i  ^11^}  -  |.V‘ 

*  r:i  c>'Ki*  +°<4i 


where  oC^,  st  0.090776  ii  the  root  of  Equation  (20)  whenir  =  m=l  . 
Equation  (30)  has  the  following  physical  interpretation.  As 
a -»  b  ,  the  cross-sectional  area  of  the  cylindrical  3hell  sub¬ 
tended  by  i  small  enough  angle  measured  from  the  axis,  may  be 
considered  as  part  of  the  cross-sectional  area  of  a  semi¬ 
infinite  slab  of  thickness  l(b-a)  .  There  is  a  symmetric 
temperature  distribution  in  the  slab  whose  isothermal  faces  have 
equal  temperature.  Hence  the  fluxes  at  the  surfaces  are  of 
equal  magnitude  but  of  opposite  sign,  showing  that  energy  flows 
out  from  both  faces.  The  dimensionless  temperature  at  the 
center  of  the  slab*  is  a  maximum  and  equal  to  I.I864. 

5  For  a  slab  of  thickness  2.  b  ,  Equation  (3)  is  $0.  * 
Frank-Kamenetzky  found  the  critical  temperature  distribution 

to  be 

ZJU  ^  [(t)1*  <31> 

where  ^=0.87846  ,  the  maximum  of 

$=2[c«&-'  (32> 

occurs  at  =UI44  which  is  the  root  of 

'  •*  «  r  ® 1  •  (33) 

6  is  the  slab  center  temperature. 
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Prom  Equations  (20)  and  (30),  Equation  (33)  may  be  obtained. 
Hence  the  temperature  defined  by  Equation  (30)  does  indeed 
correspond  To  the  slab  center  temperature.  The  critical  value 
of  S  ,  ,  for  the  slab  of  thicknes3  x  (b-a)  ,  can  be 

obtained  from  Equation  (I4.).  This  gives 


jllAfc  Sc 

ir=  1 

»Vl  -*  1 

Therefore  one  finds 


Or2)  _  (l-mf 

b1  4 


=  0.8784. 

which  has  previously  been  found  using  Equation  (3 2). 


SOLID  CYLINDER  LIMIT 


The  solution  for  the  solid  cylinder  can  be  shown  to  be  a 
special  case  of  the  more  general  hollow  cylinder  solution. 
This  is  obtained  by  setting  D=0  in  Equation  (10)  since  the 
temperature  gradient  must  vanish  on  the  axis  because  of 
symmetry  considerat ions .  Therefore  ■&  ^  o  at  ms  O  .  From 

Equation  (16)  m 


yc  =  2  ;  A  =  ■—  ( 

so  that  from  Equation  (lii). 


=  l*v4  a(-m„=JU4  c 

These  results  for  S’  ,  ,  and  »  are  indentical  with 

ChambrS's.  c  C,BUX 
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ENGINEERING  APPLICATIONS 


Of  late  there  has  been  a  trend  toward  the  design  of  large 
rocket  grains.  This  trend  toward  larger  types  has  made  it 
necessary  to  determine  how  large  a  grain  can  be  made  before  it 
presents  a  spontaneous  ignition  hazard.  It  is  unfortunate 
that  all  the  necessary  physical  constants  needed  for  this 
determination  are  available  for  only  a  very  few  propellants 
or  explosives. 

A  sample  calculation  of  the  critical  size  of  cordite 
propellant  for  surface  temperature  at  300°K  and  339°K  (150°F) 
is  given.  The  numerical  constants  for  cordite (37  are  as 
follows  r 

E=  50,000  cal  mole"1,  Z=  ID21*8  sec"1,  Q.  =  770  cal  cm" 3, 
and  A  =  5-3  x  10"9-  cal  cm"1  sec"1  deg"1.  Table  2  gives  the 
critical  radius  b  (in  meters)  which  corresponds  to  these  two 
temperatures  for  various  values  of  m  (inner/outer  diameter). 
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